Hydrogen storage of calcium atoms adsorbed on graphene: First-principles plane wave 

calculations 
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Based on the first-principles plane wave calculations, we showed that Ca adsorbed on graphene 
can serve as a high-capacity hydrogen storage medium, which can be recycled by operations at 
room temperature. Ca is chemisorbed by donating part of its 4s-charge to the empty 7r*-band of 
graphene. At the end adsorbed Ca atom becomes positively charged and the semi-metallic graphene 
change into a metallic state. While each of adsorbed Ca atoms forming the (4x4) pattern on the 
graphene can absorb up to five H2 molecules, hydrogen storage capacity can be increased to 8.4 
wt % by adsorbing Ca to both sides of graphene and by increasing the coverage to form the (2x2) 
pattern. Clustering of Ca atoms is hindered by the repulsive Coulomb interaction between charged 
Ca atoms. 

PACS numbers: 31.15.ae, 68.43.-h, 82.30.Fi 
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In order to develop an efficient medium of hydrogen 
storage, carbon based nanostructures functionalized by 
transition metal atoms have been a subject of active 
studj«!»2i2i4JL. Recently, Yoon et al& have demonstrated 
that covering the surface of Ceo with 32 Ca atoms can 
store 8.4 wt % hydrogen. Their result, which is crucial for 
safe and efficient hydrogen storage^, inspired us to con- 
sider graphene as the substrate material for Ca atoms. 
Graphene is precursor to Cgo and carbon nanotubes, but 
being a single atomic plane of graphite its both sides may 
be suitable for the adsorption of Ca atoms. Graphene 
by itself has been synthesized showing unusual electronic 
and magnetic properties^. 

In this paper, we showed that Ca atoms, in fact, can 
be bound on both sides of graphene plane and each Ca 
atom absorbing four H2 results in a medium of high- 
capacity hydrogen storage of 8.4 wt %. In the present 
case the binding energy of the fourth H2 absorbed by 
Ca atom is still significant and is ~300 meV. While each 
Ca atom donates part of its charge to the graphene layer, 
graphene, by itself, having Fermi surface consisting of six 
points at the corners of the hexagonal Brilloiun zone, is 
metallized. These results are obtained from our study 
based on first-principles calculations^. 

We first consider the adsorption of a single Ca on the 
graphene as the substrate material. This is modeled by 
one Ca atom adsorbed on the hollow site (namely HI 
site above the center of hexagon) for each (4x4) cell 
of graphene (namely one Ca atom for every 32 carbon 
atoms). The Ca-Ca interaction is indeed negligible ow- 
ing to large distance of ~ 9.84 A between them. A 
chemical bonding occurs between Ca and C atoms with 
a binding energy of 0.99 eV and Ca+graphene distance 
of 2.10 A. Similar to the bonding mechanism of Ca on 
Ceo, Ca atom donates part of its charge from 4s-orbital 
to the 7r*-bands of graphene. Due to the formation of an 
electric field between Ca atom and the graphene layer, 
part of this charge is then back-donated^ to the unoccu- 
pied 3d-orbitals of Ca through their hybridization with 



7r*-states. The resulting positive charge of Ca atom is 
calculated to be ~ 0.96 electrons^. The diffusion of the 
single Ca atom adsorbed on the graphene has to overcome 
relatively small energy barriers of Q =118 meV and 126 
meV to diffuse to the top site (i.e. on top of C atom) 
and bridge site (on top of the C-C bond), respectively. 
Ca atom adsorbed on the top or bridge sites becomes less 
positively charged (~ 0.89 and ~ 0.92 electrons, respec- 
tively) . 

A denser Ca coverage, which is energetically more fa- 
vorable, is attained, if one Ca is adsorbed on each (2x2) 
cell of graphene with Ca-Ca distance of 4.92 A. Ca atom 
adsorbed on the top and bridge sites has a binding energy 
of 0.86 and 0.89 eV, respectively. However, energetically 
most favorable adsorption site is found to be the HI site, 
which is 2.15 A above the graphene with a binding en- 
ergy of 1.14 eV. Here, the Ca-Ca coupling is subtracted 
from the calculated binding energy. In this dense (2x2) 
coverage, a stronger electric field is induced between Ca 
atoms and the graphene layer, which, in turn, leads to a 
larger back-donation of charge from the graphene layer 
to 3c?-orbitals of Ca atom. Hence by increasing Ca cover- 
age from (4x4) to (2x2), adsorbed Ca atoms become less 
positively charged, but their binding energy increases. As 
demonstrated in FiglU even if it is energetically more fa- 
vorable, the clustering of adsorbed Ca atoms is hindered 
by the Coulomb repulsion. 



We next consider the double sided adsorption of Ca. 
The binding energy of second Ca atom for the double 
sided adsorption with H1+H2 and H1+H3 configura- 
tion indicated in Fig[2] (c), is 1.27 and 1.26 eV, respec- 
tively. Since the repulsive Coulomb interaction between 
Ca atoms on the upper and lower part of the plane is 
screened by the negative charge around graphene, the 
binding energy of Ca atom in the double sided adsorption 
is larger than that in the single sided adsorption. It is 
also found that 3<i-orbitals of both Ca atoms have higher 
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FIG. 1: (Color online) (a) Top-right panel: A (4x4) cell of 
graphene having four Ca atoms. As Ca at the initial position 
is moved in the direction of the arrow, its z-coordinate is 
optimized. The remaining three Ca atoms are fully relaxed. 
Beyond the position 2 of the first Ca, the Coulomb repulsion 
pushes the second Ca atom in the same direction through 
positions 3 4 and 5 to maintain a distance with the first 
Ca. Top-left: The variation of energy as the first Ca is moves 
through positions 1-5. (b) Bottom-right panel: Two Ca atoms 
adsorbed on each (4x4) cell of graphene with their initial 
positions and . As the first Ca moves from to 1, the 
second one moves from to 1 having the Ca-Ca distance of 
3.74 A, whereby the energy is lowered by ~0.176 meV. Two 
Ca atoms are prevented from being closer to each other and as 

the first Ca moves from 1 to 2,3,4 and 5 positions, the second 

iii i 

one reverses his direction and moves through 2,3,4 and 5 
in the same direction as the first Ca atom. Bottom-left: The 
variation of the energy with the positions of Ca atoms. 



occupancies as compared with Ca atom in the single sided 
adsorption. It is noted, however, that the partial occu- 
pancy of 3d-orbitals of Ca atom does not cause any mag- 
netic properties in the system. Our results indicate that 
a stable and uniform Ca coverage up to 9=12.5 % (9= 
25 %) can be attained for single sided (double sided with 
H1+H2 or m+H3) adsorption forming a (2x2) pattern. 

Finite-temperature ab initio molecular dynamics sim- 
ulations have also been carried out for Ca adsorbed on 
the (2x2) graphene unit cell for HI geometry. Simula- 
tions are performed by normalizing the velocities of the 
ions and increasing the temperature of the system grad- 
ually from K to 900 K in 300 time steps. The duration 
of time steps are intentionally taken as 3 fs, which is 
relatively longer for a MD calculation. If the system is 
unstable, the geometry of the structure can be destroyed 
much easier in long time steps. While the bonding be- 
tween adsorbed Ca atom and the graphene layer is sus- 
tained, the adsorbed (2x2) Ca layer begins to diffuse on 
the graphene layer as the temperature of the system rises 
over ~ 300 K. However, no structural deformation is ob- 
served indicating that the Ca+Graphene system is found 
to be stable up to 900 K within 300 time steps. 

Other alkaline-earth metals, such as beryllium and 



magnesium do not form strong bonds with graphene. 
Since Be has ionization potential of 9.32 eV— which is 
much higher than that of Ca atom (6.11 eV), the charge 
of its 2s-orbital cannot easily transferred to the graphene 
layer. A similar situation occurs also with Mg having an 
ionization potential of 7.64 eV. Besides, the hybridization 
of 7T* — orbitals of graphene with the d- orbital of Ca atom, 
which is absent in both Be and Mg plays an essential role 
in strong binding of Ca to graphene. However, Ti and Co 
form strong bonds (with binding energies 1.58 and 1.20 
eV for the (2x2) adsorption pattern, respectively)^. The 
binding energies of Fe, Cr and Mo are rather weak. 



The above arguments related with the binding of Ca to 
graphene are confirmed by examining the band structure 
and the charge difference isosurfaces presented in Fig|2] 
Both H1+H2 and H1+H3 adsorption configurations are 
included in our calculations because there is a small en- 
ergy difference (H1+H2 structure is 26 meV more ener- 
getic.) between them. Hence, both adsorption configu- 
rations should be observable at room temperature condi- 
tions. Charge difference isosurfaces are obtained by sub- 
tracting charge densities of Ca and bare graphene from 
that of Ca+graphene, namely Ap — pca+Gr — Pea — PGr- 
It is seen that there is a significant charge accumula- 
tion between the adsorbed Ca atom and graphene, which 
forms the ligand field. Partial occupation of 3c?-orbitals 
of Ca can be most clearly demonstrated by the projected 
density of states in FigJ5] (b). The empty ir*— bands be- 
come occupied through charge transfer from 4s-orbitals 
of adsorbed Ca and eventually get distorted due to 3d— n* 
hybridization between 3d-orbitals of Ca and the states of 
7T*— bands as a result of the charge back-donation pro- 
cess. Occupation of distorted graphene ir*— bands gives 
rise to the metallization of semi-metallic graphene sheets 
for all adsorption sites. It is also seen that charge den- 
sity around graphene layer increased significantly as a 
result of double sided adsorption of Ca. The increase 
of charge back-donation to 3d-orbitals becomes clear by 
the increased 3d-projected density of states below the 
Fermi level. Changing the adsorption configuration from 
H1+H2 to H1+H3 does not make any essential changes 
in the electronic structure. One notes that the position 
of Fermi energy and hence electron density can be mon- 
itored by the controlled doping of Ca atoms. Metalliza- 
tion process is also important for graphene nanoribbons, 
which form conductive interconnects and spintronic de- 
vices in the same nanostructur o 12 ! 13 . It might be an in- 
teresting study to investigate the magnetic and electronic 
properties of Ca adsorption on graphene nanoribbons due 
to its different bonding mechanism. 



We next study the absorption of hydrogen molecules 
by Ca atoms. A summary of energetics and geome- 
try related with the absorption of molecular H2 for HI, 
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FIG. 2: (Color online) (a) The (2x2) cell of graphene lattice and the energy band structure of bare graphene folded to the 
(2x2) cell, (b) Single Ca atom is adsorbed on the HI adsorption site of the (2x2) cell of graphene, energy band structure 
and corresponding total density of states (dotted blue-dark curve) and partial density of states projected to Ca,-3d orbitals 
(green-gray). Isosurfaces of the difference charge density, Ap, with pink (light) and blue (dark) isosurfaces indicating charge 
accumulation and charge depletion regions. Isosurface charge density is taken to be 0.0038 electrons/A 3 , (c) Similar to (b) 
(excluding the partial and total density of states), but Ca atoms are adsorbed on both sides of graphene at the HI and H2 sites. 
(H1+H3 configuration is also shown.) (d) Partial densities of states on H-s (red-dark) and Ca-3d (green-gray) orbitals for 2, 3, 
and 4 H2 absorbed in HI configuration, and also isosurface of difference charge densities corresponding to 4H2+Ca+Graphene 
configuration. Zero of band energy is set to the Fermi energy, Ej?. 



H1+H2, and H1+H3 sites for the (2x2) and HI site for 
the (4x4) coverage are given in Figj3] The binding mech- 
anism of H2 invokes not only the adsorbed Ca atom, but 
also the graphene layer. In the case of single and dou- 
ble H2 absorption, the absorbed molecules are parallel 
to graphene and all hydrogen atoms are equidistant from 
Ca atom. As a result, both hydrogen atom of each ab- 
sorbed H2 have the same excess charge of ~ 0.08 elec- 
trons. Once the number of H2 absorbed by each Ca 
atom exceeded two, absorbed H2 molecules tent to tilt 
towards Ca atom because of increased positive charge of 
Ca atom and the symmetry of the bonding configura- 
tion of H2 molecules. The charge of Ca, H atom closer 
to Ca, H atom farther from Ca and graphene are calcu- 
lated for 8H 2 +2Ca-|-Graphene system corresponding to 

H1+H2 configuration in FigOJto be ~ +1.29, 0.06, 

~ —0.11 and ~ —1.23 electrons. One hydrogen atom of 
tilted H2, which is closer to Ca has more excess charge 
than the other one. It is important to note that charges 
transferred to absorbed H2 are not only from Ca atom. 
Graphene atoms at close proximity also supply charge 
through back-donation process. At the end, ionic bond- 
ing through attractive Coulomb interaction between pos- 
itively charged Ca and negatively charged H and weak 
van der Waals interaction are responsible for the forma- 
tion of mixed bonding between H2 molecules and Ca ad- 
sorbed on graphene. The above discussion is substan- 
tiated by the partial density of states in Fig|3fd). The 
excess charge on H-s and Ca-5d orbitals and their contri- 



bution to the states below the Fermi level increase with 
increasing number of H2 molecules. Broadening of the 
molecular level of H2 at ~ —9 eV indicates significant H2- 
H2 interaction, that in turn increases the binding energy. 
In fact, the binding energy of the first H2 molecule to Ca 
atom which prefers to be parallel to the graphene layer is 
generally small. Whereas the average binding energy for 
two H2 molecules which are again located parallel to the 
graphene layer, and for three or more H2 molecules which 
are tilted around Ca atom are larger. We note that the 
adsorption of Ca atoms and also H 2 molecules slightly 
affect the underlying graphene lattice and C-C distance. 
The average C-C distance of bare graphene is increased 
from 1.41 A upon adsorption of Ca and absorption of 
H2 to d c values indicated in all (2x2) structures in FigJ3] 
Since Ca-Ca interaction is negligible in (4x4) structures, 
there is no variation in average C-C distance. 

Maximum number of absorbed H2 per adsorbed Ca 
atom is four for the (2x2) coverage yielding a H2 storage 
capacity of 8.4 wt % and five for the (4x4) coverage of 
graphene. The reason why we include the (4x4) cover- 
age even if the resulting gravimetric density is very low 
(~ %2.3 wt), is to mimic the Ca-H.2 interaction in the 
absence of H2-H2 interaction occurring in the (2x2) cov- 
erage. Fifth H2 molecule can be bound to the top of 
Ca atom in the (4x4) coverage with a significantly high 
binding energy. Other 4 H2 molecules remain in quadri- 
lateral positions around Ca. When we compare graphene 
with C6o^, we can conclude that Cgo with a single Ca ad- 
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FIG. 3: (Color online) Sites and energetics of Ca adsorbed on graphene with the (2x2) coverage and absorption of H2 molecules 
by Ca atoms. d c is the average C-C distance in the graphene layer. El is the binding energy of Ca atom adsorbed on Hl-site, 
which is a minimum energy site. For H1+H2 or H1+H3 configurations corresponding to double sided adsorption, El is the 
binding energy of the second Ca atom and El is the average binding energy. For HI, H1+H2 and H1+H3 configurations, Ei 
is the binding energy of the first H2 absorbed by each Ca atom; E n (n=2-5) is the binding energy of the last n th H2 molecule 
absorbed by each Ca atom; E n is the average binding energy of n H2 molecules absorbed by a Ca atom. Last row indicates 
the sites and energetics of one Ca atom adsorbed on each (4x4) cell of graphene and absorption of H2 molecules by each Ca 
atom. Only the (4x4) coverage can absorb 5 H2 molecules. The shaded panel indicates energetically the most favorable H2 
absorption configuration. 



sorbed on the surface yields similar results with the (4x4) 
coverage on graphene. However, increasing Ca coverage 
adsorption results in lower binding energies of absorbed 
H 2 molecules in the present case. Unfortunately, we can- 
not comment on the case of high Ca coverage of Cgo, 
since Yoon et al& did not give details on the energetics 
of H2 molecules in denser Ca adsorption. They have just 
emphasized that adsorption of 32 Ca results in full cov- 
erage of Ceo surface and this structure can absorb up to 
92 H2 molecules with binding energy of ~ 0.4 eV. Un- 
der these circumstances, single Ca atom can hold only 3 
H2 molecules. In graphene structures, while the charge 
on (ie. charge depletion or positive charge) Ca increases 
with increased number of the absorbed H 2 molecules, the 
electric field around Ca increases. This, in turn, results in 
a decrement in the distance between adsorbed Ca and po- 
larized H2 molecules. The charge on graphene decreases, 
as well. 

In conclusion, this paper deals with two different sub- 
jects which are of current interest; namely graphene and 
hydrogen storage. Firstly, we showed that recently syn- 
thesized graphene with carriers behaving as if massless 
Dirac fermions can be metallized as a result of the ad- 



sorption of Ca atoms. Electrons donated by Ca is ac- 
commodated partly by the 7r*— bands of graphene, partly 
back donated to its 3d-orbitals. Ca atoms can be bound 
to both sides of graphene and can attain 25% coverage 
without clustering. Secondly, we found that each ad- 
sorbed Ca can absorb up to four hydrogen molecules. At 
full coverage this yields a storage capacity of ~8.4 wt%, 
which is higher than the value set for feasible gravimet- 
ric density of hydrogen storage. The calculated bonding 
energies of hydrogen molecules are suitable for room tem- 
perature storage; while above the room temperature hy- 
drogen molecules are released, Ca atoms remain adsorbed 
on graphene for further recycling. Even though storage 
capacities higher than the present case is achieved^ in 
different nanostructures, our results may be important 
for efficient hydrogen storage since graphene flakes are 
now easily available. 

This work is supported by TUBITAK through the 
grant TBAG104536. Part of the computational resources 
for this study has been provided through Grant No. 2- 
024-2007 by the National Centre for High Performance 
Computing of Turkey, Istanbul Technical University. 



5 



* Electronic address: ciraci@fen.bilkent.edu.tr 

1 S. Dag, Y. Ozturk, S. Ciraci and T. Yildirim, Phys. Rev. 
B, 72, 155404 (2005). 

2 Y. Zhao, Y.-H. Kim, A.C. Dillon, M.L. Heben, and S.B. 
Zhang, Phys. Rev. Lett. 94, 155404 (2005). 

3 T. Yildirim, S. Ciraci, Phys. Rev. Lett. 94, 175501 (2005). 

4 E. Durgun, S. Ciraci, W. Zhou, T. Yildirim, Phys. Rev. 
Lett. 97, 226102 (2006). 

5 C. Ataca, E. Akturk, S. Ciraci, H. Ustunel, Appl. Phys. 
Lett. 93, 043123 (2008). 

6 M. Yoon, S. Yang, C. Hicke, E. Wang, D. Geohegan, Z. 
Zhang, Phys. Rev. Lett. 100, 206806 (2008). 

7 R. Coontz and B. Hanson, Science 305, 957 (2004). 

8 K. S. Novosolov and A.K. Geim, Nature 438, 197 (2005). 

9 We have performed first-principles plane wave calculations 
within Density-Functional Theory using PAW potentials 
[P.E. Blochl, Phys. Rev. B 50, 17953 (1994)]. We used 
Local Density Approximation (LDA), since the van der 
Waals contributions are generally better accounted. We 
also calculated binding energies by using Generalized Gra- 
dient Approximation (GGA). We found that GGA binding 
energy for n = 1 is 120 meV smaller than the LDA binding 
energy for HI adsorption site. However for n = 2, the bind- 
ing energy becomes ~ 210 meV smaller than LDA because 



the H2-H2 interaction is neglected in the former. Hence 
GGA and LDA result in very close atom configurations. 
Numerical results are acquired by using VASP [G. Kresse 
and J. Hafner, Phys. Rev. B 47, R558 (1993)]. A plane- 
wave basis set with kinetic energy cutoff ft 2 |k + G| 2 /2m 
= 900 eV has been used. The Brillouin zone has been 
sampled by (19x19x1) and (9x9x1) special mesh points in 
k-space for (2x2) and (4x4) graphene cells, respectively. 
In the course of structure optimizations, the convergence 
for energy is chosen as 10 -6 eV between consecutive steps, 
and the maximum force allowed on each atom is less than 
10" 2 eV/A. The pressure on the system is kept smaller 
than ~ 0.1 KBar per unitcell. 

Charge transfer values are calculated according to the 
Bader analysis. See G. Henkelman, A. Arnaldsson, H. Jon- 
sson, Comput. Mater. Sci. 36, 254-360 (2006). 

1 C. Kittel, Introduction to Solid State Physics, 7th ed. (Wil- 
ley, New York, 1996). 

2 H. Sevincli, M. Topsakal, E. Durgun and S. Ciraci, Phys. 
Rev. B, 77, 195434 (2008). 

3 Y.-W. Son, M.L. Cohen and S.G. Louie, Nature, 444, 347 
(2007); H. Sevingli, M. Topsakal and S. Ciraci, Pyhs. Rev. 
B, 78, 245402 (2008). 



